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Abstract 

In this paper we consider asymptotic behavior of a hybrid action of f(R) gravity 
model which proposed by Saffari & Rahvar (2008), in the Solar system scale, which 
can explain Pioneer anomalous acceleration. We use the resultant weak field grav- 
itational potential which comes from the hybrid action to test its impacts on the 
Solar system dynamics, by comparing theoretical precession of perihelion of a test 
particle, cb, with corrections to the standard Newtonian - Einstenian precessions 
of perihelia of some planets, which recently estimated by Pitjeva (2005a,b, 2006, 
2008). Here we show that the asymptotic behavior of hybrid action, are in more 
accordance with observation relative to the other modifications such as power law 
and logarithmic corrections (Iorio, 2008). We also show that an extra additional 
lensing of the prediction of General Relativity is reproduced. Finally we consider 
the stability condition of planetary orbits in presence of the hybrid action. 
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1 Introduction 



To accommodate some recently observed phenomena (jRiess et all 11998 



2006 



Perlmutterl.119991 IColeLl2005l iTegmarkL l200d ISpergelL 120031 . 120071 : Idowe et al 



Anderson et all Il998l 120021 ). occurring at very different scales ranging 



from Solar system to cosmological distances, which at present, have not yet 



* Corresponding author 

Email address: rsk@guilan.ac.ir (Reza Saffari). 



Preprint submitted to Elsevier 



14 December 2010 



found fully satis factorily explana t ions in terms of conven t ional physics, grav- 



itation or not, ( Bertolami et al.l . [2006: B ertolamil . l2008t iLammerzahl et al. 



20081 ). we considering modifications of the standard laws of gravity at different 



scales; As well as done in recent years. These anomalous effects are: The recent 



data coming from the luminosit y distance of SuperNovae la (SNela) (jRiess et al. 



19981 ). wide galaxy surveys (IColel. 120051 ) and the a nisotropy of cosmic mi- 



crowave background radiation ( jSpergell . 120031 120071 ) suggest that about 70 
percent of the energy density of the present universe is composed by dark en- 
ergy, responsible for an accelerating expansion. Further more in small scales, 



the d ata coming from galactic rotation curves of spiral galaxies (I Clowe et al 



20061 ) can not be explained on the basis of Newtonian gravity or General 



Relativity (GR), if one does not introduce invisible dark matter. In Solar sys- 
tem there is an unexplained acceleration a p = (8.47 ± 1.33) x 10~ 10 ms -2 ap- 
proximately directed towards the S un affecting the Pioneer 10/1 1 spacecrafts 
after they passed 20 AU threshold (j Anderson et all 1 19981 . 120021 ). In order to 
reconcile theoretical models with observ ations, existence of dark matter is 
postulated in Solar system (jNietol . 120081 ) . The main problem of dark energy 
and dark matter is understanding their nature, since they are introduced as 
ad hoc gravity sources in a well define model of gravity. Therefore we find 
another possibilities, like various modification of gravity. In this paper we 
are concerned with f(R) gravity theory, where th e gravitational lagr a ngian 
depends on a function of t he scalar curvature R (INojiri fc Odintsovl . 120071 : 



Capozziello fc Francaviglial . l2008t iFaraonil . 120081 ). These theories are also re- 



ferred to as 'extended theories of gravity', since they naturally generalize GR: 
in fact, when f(R) = R the action reduces to the usual Einstein-Hilbert action, 
and Einstein's theory is obtained. These theories can be studied in the metric 
formalism, where action is varied with respect to metric tensor and in Palatini 
formalism, where action is varied with respect to metric and affine connec- 
tions, which are supposed to be independent from one to another. In general, 
two approaches are not equivalent: th e solutions of P alatini field equations 
are a subset on metric field equations (iMagnanol . ll994T ). Actually, f(R) theo- 
ries provide cosmologically viable mode ls, where both the inflation phase and 
the accelerating expansio n are reported ( de la Cruz-Dombriz fc Dobadol . 2006 : 



Nojiri fc Odintsovl . l2009l ). Furthermore, they have been used to explain the 



rotat ion curves of galaxies without need for dark matter (IFrigerio fc Saluccil . 



20071 ) . However, because of the excellent agreement of GR with Solar system 



and binary pulsar observations, every theory that aims at explaining galaxies 
dynamics and the accelerating expansion of the Universe, should reproduce 
GR at the Solar system scale, in a suitable weak field limit. Without going 
into the details, we want to test the model of f(R) gravity and using perturba- 
tive approach, we focus on the impacts which modification of the gravitational 
field due to nonlinearity of the gravitational lagrangian, have on Solar system 
dynamics. We compare the theoretical prediction for the precession of the lon- 
gitude of the pericenture of a test particle with the correction to the standard 
Newtonian-Einsteinian precession of perihelia of some planets recently esti- 
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mated by Pitjeva.t lPitjeval . l2005al Jbl). At first, we derive an explicit expression 
of the secular, averaged over one orbital revolution perihelion precession, in- 
duced by anomalous acceleration of asymptotic behavior of hybrid action on 
the orbits of Solar system planets. Then we will compare obtained results in 
small eccentricity approximation, with the latest estimated corrections of the 
perihelion rates. It is important to note t hat the corrections to the perihelion 
rates determined in (jPitjeval . 120061 . 120081 ). are phenomenologically estimated 
quantities of a global, least-square solution in which only Newtonian and Ein- 
stenian dynamics was modeled; no exotic dynamical terms were included in 
the fit. Thus in our opinion such phenomeno logical corrections can genuinely 
be used to get information on a hypothetic, un-modeled force. Also, central to 
the development of the good modified theory of gravity that can describe the 
phenomenology usually ascribed to dark matter in the analysis of the amount 
of deflection of light implied by the theory. Then we obtain deflection of light 
for the hybrid model and compare it with the prediction of GR. This paper is 
organized as follows: in section [21 we briefly review the theoretical formalism 
of f(R) gravity, in section [3j we outline a general approach to the perturba- 
tions of the gravitational field of GR, in section HI we introduce asymptotic 
behavior of hybrid f(R) model and we compare theoretical predictions with 
observations and compare our results with some models that tested in Solar 
system, in section [5j we calculate deflection of light, in section [6] we consider 
stability of circular orbits and finally, conclusions are in section [7J 



2 f(R) gravity 



He re we introduce the field equation of f(R) gravity (ISotiriou fc Faraoni 
20101 ). we start from modified Einstein-Hilbert action 



S = -r 

2k 



J ^x^f(R) + S n 



The gravitational part of the Lagrangian, f(R) is represented by a function 
of the scalar curvature, R. S m is the matter part of the action. In metric 
formalism affine connections, T's are supposed to be Levi-Civita connections 
of g and consequently, the scalar curvature R has to be intended as 

R = R{g)=gf*Rap(g). (2) 



On the contrary in the Palatini formalism metric g and affine connections are 
supposed to be independent, so that the scalar curvature R has to be intended 

as 

R = R{T)=g^R a p(V), (3) 
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where R Q p(r) is the Ricci-like tensor of the connections T. We follow the 
metric formalism that in this regime, Eq. ([I]), is varied with respect to the 
metric g^, and obtain the field equation of motion 

F(R)R^ - V M V,F(i?) - kT^ = g^f(R) - V Q V a F(i?)), (4) 



where is the standard minimally coupled matter energy- momentum tensor 
and F = df jdR. Contraction of the field equation Eq. (j4]), with metric tensor 
leads to scalar equation 

3V Q V a F(R) + F(R)R - 2f(R) = ^-T, (5) 



where T is the trace of the energy-momentum tensor. Eq. ([5]) is a differential 
equation for the scalar curvature R. In order to compare the prediction of 
f(R) gravity with Solar system data, we have to consider the solutions of the 
field equation Eq. (T4j) supplemented by the constraints Eq. ([5]) in vacuum, 
since tests are based on the observations of the dynamics of the planets in the 
gravitational field of the Sun. We notice that if R = constant we obtain the 
Palatini case: so for a given f(R) function in vacuum, solutions of the field 
equation of P alatini f(R) gravity are a subset of the field equation of metric 
f (R) gravity. flMultamaki & Vilial . l2006h 



3 The proposal method to test modified f(R) gravity model 

In this section we outline the general procedure that we are going to apply 
to the solution of f(R) gravity, In order to compare the prediction of these 
gravity models with the existing data. In weak field approximation of metric 
and spherically symmetrical condition 

ds 2 = -B(r)dt 2 + A{r)dr 2 + r 2 sin 2 8d<p 2 + r 2 d8 2 , (6) 

the gravitational (scalar) potential 0(r) is read from the B(r) function 

B(r)-1 + S£i. (7) 

We expect a gravitational potential in the form of 

0(r) = <p N {r) + A0(r), (8) 
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-GM/r is the Newtonian potential of a point-like mass M in 

lN i 



where (f) (r) 

center, and A0(r) <C (r) is a correction vanishing for /(i?) — )■ -R. Now if 
we have a correction for Newtonian potential, we can calculate the perturb- 
ing acceleration and then we will be able to cal culate its ef fect on planetary 
motions, within standard perturbative schemes ((RoyJ, |2005|) . Gauss equation 
for entirely radial perturbing acceleration A r is 



dzu 
~dt 



nae 



-A r costp. 



(9) 



After being evaluated onto the unperturbed Keplerian ellipse, the acceleration 
A r must be inserted into Eq. ([9]); then, the average over one orbital period 
P must be performed. To this end, we use following equations; where n = 
GM/a 3 , is the Keplerian mean motion, a is the planet s semi major axis, e 



is the eccentricity and ip is the true anomaly. (jRoyj, 120051 ) 



COS if 



cos E 



1 — e cos E ' 



(10) 



a(l — ecos-E), 



1 — e cos E , ^ 
dt = dE. 

n 



(11) 
(12) 



In fact, what we aim at, is evaluating the perturbations induced on longitudes 
of perihelia by the corrections to the gravitational field due to f(R) gravity. 
In orde r to com pare them with the latest observation we use table 1 and 2 of 
Pitjeval ( l2005al lbl). that recently processed almost one century of data of dif- 
ferent types for major bodies of the Solar System in the effort of continuousl y 
improving the EPM2004:/ EPM200Q planetary ephemerides, (jPitjeval 120061 ). 
Among other things, she also estimated corrections to the secular rates of the 
longitudes of perihelia of some planets of the Solar system as fit-far parameter 
of a global solution in which she contrasted, in a least-square way, the obser- 
vations to their predicted values computed with a complete set of dynamical 
force models including all the known Newtonian and Eisteinian features of 
motion. As a consequence, any un-modeled exotic force present in nature is; 
in principle, entirely accounted for by the obtained apsidal extra precession 
Azu. See table [TJfor inner planets and table [2]for outer planets of Solar system. 



Now let us briefly outline how we are going to put f(R) gr avity on the test. 
We follow Iorio's regime for testing modified gravity models, (jlorioj, 120081 ) and 
take the ratios of observational Aw for different pairs of planets A and B 



n 



AB 



Am a 
Azupt ' 



(13) 
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and compare it with the predicted ratios of modified model of gravity £ab- If 
we have ^ ab — \^ab — £ab\ — within the errors, the f(R) gravity model ex- 
amined, can still be considered compatible with data, otherwise it is seriously 
challenged. 



4 Asymptotic behavior of the hybrid action 



Saffari fc Rahvarl ( 120081 ). use the inverse solution to extract an appropriate 



action for modified gravity, 

R + A 



f(R) = R + A + 



In 



R + A 



R/Ro + 2/a R, 



(14) 



where Rq = 6a 2 /d 2 and R c is the constant of integration, a ~ 10~ 6 , is a small 
dimensionless constant and d ~ lOkpc, is length scale in the order of galactic 
size and A is cosmo logical constant. In Solar system, for the range of R ^> A 
and R/Rq 3> 2/a, action reduces to 

f{R) = R + R In (15) 



On the other hand, in galactic scales, for cr C 1 and R ~ Rq ~ A generic 
action can be written as 

/(i?) = ( J R + A)[l + |ln(^)], (16) 



which can justifies flat rotation curves of spiral galaxies, and for small a can 
be written as 

(R + A) 1+a / 2 

tic 



which reduces to f(R) = R + A for a <C 1. In Solar system scale, in which the 
range of distances are sufficiently smaller than the length scale of the hybrid 
action, we obtain metric elements of Eq. (jSj) up to the first order terms in a 
as 



. 1 2GM ar 
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therefor the weak field effective potential which comes from Eq. (fl5j) obtains 

r, (19) 



as 

GM ac 2 



2d 



where the first term is the Newtonian potential and the second one is corrected 
term. Acceleration of a test particle in this potential is 

GM «" (20) 



r 2 2d ' 

where the second term in the right hand side of this equation, is a constant 
acceleration, independent of the central mass. We may correspond this extra 
term to the pioneer anomalous and constrain it with the observed value of 
a p = (8.47 ± 1.33) x l(T 10 ms _2 which results in a/d = 10-26m -1 . In this 
case, correction to the gravitational potential is 

A0(r) = ^r, (21) 



where the acceleration of a test particle in this potential is 

A - -°4. (22) 



On using 



2t 



J (cos E - e)dE = -2vre, (23) 



and Eqs. (191- (1T21) it yields to the following perihelioin precession 



<zu >= — ; (24) 

2d y/GM V ; 



It is important to note that w depends on a positive power of the semi major 
axis a. The predicted extra precession of Eq. (1241 can be fruitfully compared to 
the correction of usual Newtonian-Einsteinian p erihelion rates o f t he planets o f 
the Solar system phenomeno logical estimated by lPitjeva d2005al fbL I200H . l2008h . 



According to the general outline of section EJ we will not use one perihelion at 
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a time for each planet. Indeed let us consider a pair of planets A and B, take 
the ratio of prediction of Eq. (12"4"|) 



W B 



1 - e A )M B a A 
\{l-e B )M A a B 



(25) 



Now we consider a pair of planets A and B and take the ratio of their esti- 
mated extra-rates of perihelia in approximation small eccentricity; if Eq. (1241) 
is responsible for them, 

,Tr i aA \-\ ?nr\ 

^ab = \- ( — ) 2 |> (26) 

w B a B 



then the quantity must be compatible with zero. Table [3] shows the results 
of observational ratio and theoretical ratio of planets. It is important to note 
that the uncertainty in 9 ab has been conservatively estimated by linearity 
adding the individual terms coming from the propagation of the errors in w 
and a in Eq. ( l26l) . In general, deviation of each function has the form of 

Qt = /(EW(a) J . ( 27 ) 

and errors in Qab due to Sa and 8 Aw have the form below 

69 ab < |^|(^4 + + + — )■ (28) 

AziiD Arij/i Atijn 2 an cla &r 



To see that how, this model acts in the Solar system, we can compare it with 
other corrections in Solar system. To this end we use power - law a nd loga- 
rithmic model in weak field approximation. ICapozziello et al.l (120071 ) started 
the Lagrangian of the form f(R) = foR n , i n the metric approach, to obtain 
solutions to describe galaxies rotation curves without need for dark matter. 
They obtain modified gravitational potential of the form 



r r r 



(29) 



where (3 relates to the the power of n in the modified Lagrangian via 



12n 2 -In- + 12n 3 - 83n 2 + 50n + 1 

^~ Qn 2 -An + 2 ' ^ ' 
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This equation clearly leads to the radial acceleration 
_ (g - l)GM 

fi r — a I 



(31) 



It yields the following perihelion precession ( lloriol . 120081 ) 



< w > 



(32) 



Capozziello et al.l ( 120071 ) applied Eq. (129|) to a sample of 15 low surfaces bright- 
ness galaxies with combined HI and Ha measurements of the rotation curve 
extending in the putative dark matter dominated region. They obtained a 
good agreement between the theoretical rotation curves and the data using 
only stellar disk and interstellar gas when the slope n of the gravity Lagrangian 
is set to the value n = 3.5 (giving = 0.817) obtained by fitting the SNela 
Hubble diagram with the assumed power law f(R) model. For (3 = 0.817, 
results are in figure [TJ Logarith mic correction suggested to solving the prob- 
lem of dark m atter in galaxies (IVan Moorsell . Il987t iFabris fc Camposl . 120081 : 
Soboutil . 120071 ) . The potential is the form of 



V ln = -aGM ln(- 



r 



(33) 



that gives the radial extra acceleration 
aGM 



(34) 



and it gives the precession to the form of 




(35) 



The results for this model are in figure [TJ As a conclusion the log arithmic an d 



power law correction is ruled out by the present day observations ( lloriol . 120081 ) . 
As we see in figure [TJ the correction to the form of this asymptotic behavior 
of hybrid f{R), has the minimum difference with observational data. 

Now the question is, about the properties of these models, we see that each 
of these corrections, has the term like r n in it. For power law model n = 0.817, 
for logarithmic correction n = — 1 and for the hybrid model n = 0. then we can 
debate on the power of r n , and find that, which power of n has the appropriate 
accordance with observational data. With primary surveys we obtain n = 5.38. 
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For such potential, we can find the metric and the f(R) model from inverse 
solution which is in preparation. 



5 Deflection of light 



To see that the effect of this asymptotic hybrid f(R) model on the amount 
of deflection, we calculate the deflection of the orbit from a straight line 

A$ = 2|$(r ) - $ 00 | - 7r, 
where r is the closet distance to the Sun. 



A$ 



oo 

= 2 J A^ir) 



' B(r ) 
B{r) 



-1/2 



dr 



7T. 



By using the values of A(r) and B{r), we could obtain A$ as elliptic integrals 



A$ = GR terms - ~r \in - + In ./(-) 2 - 1] 



+i^Hr In - + In /(-)* - 1] - J(-) 2 " !]■ (36) 

We replace r with r , then to first order in MG/r the deflection is 

A$ = 37) 



and this is the prediction of GR. Then this asymptotic behavior of hybrid 
f(R) model do not violate GR regimes and have no any additional effect on 
deflection of light. 



6 Stability of circular orbits 



Finally in the context of stability condition of orbits, we test the obtained 
weak field potential with numerous equation of stability: 

irA + - > 0, 38 
F c (r) r 
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in which F c is central force, we obtain r 2 < GMd/3ac 2 . Then maximum of 
the radius of stable circular orbit, is about r s 2.22 x 10 11 km which is 10 3 
larger than typical Solar system radius, in order of magnitude. Therefore all 
the planets in this modified regime may have stable orbits. 



7 Comments and conclusions 



In this paper we have studied the secular precession of the pericenture of 
a test particle in motion around a central mass M whose Newtonian gravi- 
tational potential exhibits a correction which is to form of Eq. ( TT9|) . In order 
to put on the test the hypothesis that such extra forces are not zero we de- 
vised a suitable test by taking into account t he ratios of the c orrect i ons t o 



the secular precession of perihelia estimated by iPitjeval (l2005aU bl. 12006k 120081 ) . 
The results obtained, resumed by Eq. fT2"4"|) . show that modifications of the 
Newtonian potentials like this examined in this paper are compatible with the 
currently available apsidal extra-precession of the Solar system planets. Of 
course we have to say that this model will be ruled out if we take the inverse 
ratio of precession for planets. Hence we take this regime that, take the inner 
planet to outer planet. Our correction can be used to describe the problem of 
dark matter in the Solar system, if we choose the appropriate values for pa- 
rameters of the model. It gives Pioneer anomalous acceleration for the suitable 
value of distance. When we use the other estimative corrections to Newtonian- 
Einstenian planetary perihelion rates, maybe obtain better compatibility for 
this model in the Solar system. We had shown that this asymptotic hybrid 
model give the GR, deflection of light and had no extra effect on it. From 
comparison alternative models, we could offer the generic model of potential, 
and obtain it's action. This proposal model has minimum variance, and our 
aim is checking its compatibility with dark matter and dark energy data set. 

Acknowledgments: We would like to thank anonymous referee for useful 
comments. 
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Table 1 

observational data for inner planets 



planet Mercury Earth Mars 

A±(arcs/cy) -0.000036±.005 -0.0002±.00004 0.0001±.0005 

a{AU) 0.387 1 1.523 

e 0.205 0.016 0.0934 

P(yr) 0.24 1 1.88 



Table 2 

observational data for outer planets 



planet 


Jupiter 


Saturn 


Uranus 


Avj(arcs/cy) 


0.0062±.036 


-0.92±2.9 


0.57±13 


a{AU) 


5.203 


9.537 


19.191 


e 


0.0483 


0.0541 


0.0471 


P(yr) 


11.86 


29.45 


84.07 



Table 3 

results, A B denotes the pair of planets used £Iab 



The 

limb *i — — — 1 i' Qfl ' i 

perihelion extra rates for the planets have been retrieved from lPitieval (|2006l ). their 
errors are not formal; and re scaled by a fact or 10. The unce rtainties in the semi 
major axis have been retrieved from table 5 of iPitieva (|2005ah 



Pairs 


A 


B 




Cab 






^AB 


1 


Mer 


Jup 


-0.6±4.1 


0.272± 


10" 


-9 


0.87±4.1 


2 


Ear 


Jup 


-0.03±0.25 


0.438± 


10" 


-8 


0.46±0.25 


3 


Mar 


Jup 


0.02±0.17 


0.541± 


10" 


-9 


0.52±0.17 


4 


Mer 


Sat 


0.004±0.017 


0.201± 


10" 


-8 


0.19±0.017 


5 


Ear 


Sat 


0.0002±0.0011 


0.323± 


10" 


-10 


0.32±0.0011 


6 


Mar 


Sat 


-0.0001±0.0009 


0.399± 


io- 


-8 


0.399±0.0009 


7 


Jup 


Sat 


-0.006±0.060 


0.738± 


io- 


-8 


0.73±0.060 


8 


Mer 


Ura 


-0.006±0.152 


0.142db 


10" 


-8 


0.146±0.152 


9 


Ear 


Ura 


0.0003±0.0087 


0.228± 


10" 


-9 


0.21±0.0087 


10 


Mar 


Ura 


0.0002±0.0048 


0.281± 


io- 


-8 


0.27±0.0048 


11 


Jup 


Ura 


0.01±0.31 


0.52± 10" 9 


0.51±0.31 



14 



15r 
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Fig. 1. This figure shows accordance of observational precession with hybrid action 
theoretical results against power-low and logarithmic corrections. 



15 



